Unconventional vortex states in nanoscale superconductors due to shape-induced 
resonances in the inhomogeneous Cooper-pair condensate 
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Vortex matter in mesoscopic superconductors is known to be strongly affected by the geometry of 
the sample. Here we show that in nanoscale superconductors with coherence length comparable to 
the Fermi wavelength the shape resonances of the order parameter results in an additional contri- 
bution to the quantum topological confinement - leading to unconventional vortex configurations. 
Our Bogoliubov-de Gennes calculations in a square geometry reveal a plethora of asymmetric, giant 
multi-vortex, and vortex-antivortex structures, stable over a wide range of parameters and which 
are very different from those predicted by the Ginzburg-Landau theory. These unconventional states 
are relevant for high-T c nanograins, confined Bose-Einstein condensates, and graphene flakes with 
proximity-induced superconductivity. 
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In the last decades, the effect of the boundary on meso- 
scopic superconductors with dimensions comparable to 
the penetration depth A or the coherence length £ has 
been intensively studied [lj-Q . In applied magnetic field, 
it was found that the vortex states strongly depend on 
the size and geometry of the sample and are generally 
different from the Abrikosov triangular lattice observed 
in bulk type-II conventional superconductors (where only 
the vortex- vortex interaction plays a role). For example, 
a giant vortex induced by strong boundary confinement 
was predicted Q as the ground state in disks which was 
subsequently observed experimentally 0, [Io[ • In square 
samples, a peculiar state with an anti- vortex at the center 
surrounded by four vortices was predicted theoretically 
111 Il2| , but never observed experimentally up to now. 

All of the above theoretical works are based on the 
Ginzburg-Laudau (GL) theory. When the superconduc- 
tor is downscaled to nano-meter sizes, quantum confine- 
ment [131 ] leads to unique phenomena, especially in sam- 
ples with dimensions of the order of the Fermi wave 
length Xp. The GL theory is no longer applicable in 
this regime and the microscopic Bogoliubov-de Gennes 
formalism is required. The discretization of the energy 
levels around the Fermi level Ep was shown to induce 
quantum-size effect fl4 . 15|. quantum-size cascades (l6j 
and the shell effect [17| . As one of the important results, 
Ref. [15| reported the wave-like inhomogeneous spatial 
distribution of the order parameter, further enhanced at 
the boundary due to quantum confinement. The latter 
is important since it is well known that vortices tend to 
migrate and bepinned in areas where superconductivity 
is suppressed [l8(, i.e. it is energetically favorable for a 
vortex to suppress the superconducting order parameter 
in the region where it is already weak. In reality the be- 
havior is much more complex and in some instances the 
vortex can be pinned where the gap is large |l9j . The 
appearance of oscillations in the order parameter profile 
due to quantum confinement is thus expected to influ- 
ence the vortex states. For conventional superconduc- 



tors, kp^o ~ I0 3 {Uf is the Fermi wave vector and £o is 
the BCS coherence length), systems of size comparable 
to Xp will not be large enough to host a vortex (being 
much smaller than the coherence length). However, ma- 
terials with small coherence lengths, e.g. high-T c cuprate 
superconductors, will have kp £o wl-4 and therefore in 
such systems it is possible to obtain vortex states in the 
quantum confinement regime. Another such system is 
a graphene flake deposited on top of a superconductor. 
Because of the proximity effec t, supe rconducting corre- 
lations will diffuse in graphene [20H23l |. Such a system is 
in the clean limit since the scattering length in graphene 
is large. More importantly, in graphene, near the Dirac 
point, the Fermi wavelength is very large and can be eas- 
ily manipulated by doping. In other words, kf^o can be 
tuned, which will allow for different vortex patterns to 
be realized in the graphene flake in the quantum confine- 
ment regime, but for more accessible sample sizes (above 
lOOnm). 

In order to experimentally detect vortex states in nano- 
sized superconductors, one can extract information about 
the local density of states (LDOS) from measurements of 
the differential conductance with scanning tunneling mi- 
croscopy (STM) fiol . 24, [25 1 . An extensive analysis of the 
LDOS profile of the vortex states has been performed in 
the past [26M30| | . It is generally known that the bound 
states in the vortex core lead to peaks in the LDOS at 
energies below the superconducting gap, though the ex- 
act formation of peaks in the spectrum of a multiple flux 
line (giant vortex) will depend on the vorticity [31(. Also, 
when kp£o is small, the spectrum becomes particle-hole 
asymmetric and the lowest vortex bound state has a fi- 
nite energy [2H ]. In the quantum confinement regime, 
there exist strong vortex- vortex and vortex-boundary in- 
teractions and the quasiparticle spectrum becomes much 
more complicated. In this case the lowest bound state 
peak position does not generally coincide with the vor- 
tex core [13]. Furthermore, in case of strong interactions, 
vortex and surface bound states may combine to form a 
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complex state where LDOS contributions of individual 
constituents are not clearly visible. 

In this Letter we report novel vortex states that ap- 
pear from the interplay between quantum confinement, 
inhomogeneous superconductivity, an external magnetic 
field, and the sample geometry, in a nano-sized supercon- 
ducting square. We performed calculations for a sample 
in the quantum limit by solving Bogoliubov-de Gcnncs 
(BdG) equations self-consistently. In what follows, we 
keep constant the size of the sample and the bulk coher- 
ence length £o = hvp/nAo (where vp is the Fermi veloc- 
ity and Ao is the order parameter at zero temperature), 
while we change the parameter kp^Q and thereby tune 
the influence of the confinement on the vortex structure. 
We start from the well-known BdG equations: 

[K - E F ] u n C?) + A(l*)t; n (l*) = E n u n {l>), (1) 
A(^)* Un (^) - [K* - E F ] v n {l>) = E n v n {l>), (2) 

where Ko = (iKV + e~A*/c) 2 /2m is the kinetic energy 
and Ep is the Fermi energy, u n {v n ) are electron(hole)- 
like quasi-particle eigen-wave functions, E n are the quasi- 
particle eigen-energies, and A is the vector potential (we 
use the gauge V • A = 0). 

The pair potential is determined self-consistently from 
the eigen-wave functions and eigen-energies: 

A(?)=g «n(^)<(^)[l-2/ n ], 

E n <E c 



(3) 



where g is the coupling constant, E c is the cutoff energy, 
and /„ = [I + cxp(_E„/fcsT)] _1 is the Fermi distribution 
function, where T is the temperature. We consider the 
two-dimensional problem and assume a circular Fermi 
surface. The confinement imposes Dirichlet boundary 
conditions (i.e. u n [~?) = 0, = 0, r G dS) such 

that the order parameter vanishes at the surface. In an 
extreme type-II superconductor (and/or very thin sam- 
ple), it is reasonable to neglect the contribution of the 
supercurrent to the total magnetic field. For such a case, 
we discretize Eqs. ([TJ3]) and by using the finite difference 
method we solve them self-consistently. 

The free energy [33, 34] of the system is then calculated 
as: 

F = J2 2E nfn + k B T[f n In /„ + (f - /„) ln(I - /„)] 



- / d~V 
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J2 - 2/„) «»'<'(! - 2/nO 



(4) 



where the spatial dependence of u n and v n is implicit. 
The local density of states N(r,E) is calculated from 

N(r,E) = -Y}.f'{E n -E)\u n \ 2 + f'{E n +E)\v n \ 2 ]. (5) 
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FIG. 1. (Color online) Free energy as a function of the 
magnetic flux through the square sample, for (a) £;f£o = 2, 
and (b) fc^Co = 3. Here, Fq — h 2 /2rn^Q. The insets show 
the contour plots of the order parameter with the diagonal 
profiles in the absence of applied magnetic field. 



In this paper, we consider as an example a thin super- 
conducting square with size 5£o x 5£o- The microscopic 
parameters are set to keep Aq/E c — 0.2. The calcula- 
tions are done for different parameters kp^Q. Since we 
consider the zero temperature case, the system is always 
in the quantum limit (where T < l/kp^o). Fig. [T]shows 
our numerical results for the free energy of the found 
stable vortex configurations (the states with up to five 
vortices are shown) for two values of the key parameter, 
^fCo- The insets show the inhomogeneous profile of the 
superconducting order parameter in the absence of an 
applied magnetic field which is expected to strongly in- 
fluence the vortex structure. When comparing with con- 
ventional free energy curves obtained from the GL theory 
[35( 1 . many differences can be observed. First, the pene- 
tration field for the first vortex is suppressed because the 
order parameter is not homogeneous, allowing the vortex 
to penetrate easier at locations where the order param- 
eter is weakened. Second, the stability range in flux for 
different vortex states (with vorticity L) is not mono- 
tonically decreasing towards I$o a s L increases. More- 
over, those stability ranges strongly vary when kp^o is 
changed! For example, for fc^^o = 2 the vortex structures 
with even vorticity are stable over a broader magnetic 
field range while for kp^o — 3 surprisingly the structures 
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FIG. 2. (Color online) Transition fluxes in units of $0 between 
ground states with consecutive vorticities for different /cf£o- 



with odd vorticity are the favored ones! The main reason 
behind this phenomenon is that different confinement- 
induced oscillations in the order parameter for different 
k F S,o value will stabilize different symmetries of the vor- 
tex pattern. To illustrate this effect better, we plot in 
Fig. [5] the applied magnetic flux at which ground-state 
transitions between states with consecutive vorticities oc- 
cur, as a function of /cf?o- Notice the varying ranges of 
stability of different vortex states, which are very sen- 
sitive to k F ^Q. Of course, for large k F £o the behavior 
of the system converges to a more conventional picture 
(with each new vortex entering the system with roughly 
one flux-quantum added). 

To underpin the reasons for this varying stability of 
vortex states in what is otherwise a rather simple, square 
system, we show in Fig. [3] some of the typical states for 
the case of k F £ = 2 (the order parameter, its phase, and 
corresponding LDOS). As elaborated above, the quan- 
tum confinement of electrons here strongly affects the 
spatial distribution of the order parameter [see inset in 
Fig.QJa)], having three oscillations across the square and 
four distinct minima that enhance the fourfold symme- 
try. This automatically leads to the improved stability of 
states with even vorticity, similar to the case of a square 
with 4 antidots ;36]. We also observe that, before ceasing 
at the boundary, the order parameter is enhanced near 
the surface, with the highest value found near the cor- 
ners. Due to the effect of the boundary and the shape 
resonances, the order parameter is also enhanced at the 
center of the square . To reiterate a fairly obvious point, 
vortices are repelled by the peak positions of the super- 
conducting pair amplitude, and the four low amplitude 
locations (with value only 2/3 of the peaks) will pin vor- 
tices rather strongly. Figs.|3{a-c), show the L — 1 ground 
state for applied flux $/$o = 4. Surprisingly, we find 
that the only vortex in this state is actually sitting in 
one of the minima of the order parameter and the four- 
fold symmetry is broken. We emphasize that this state 
is not possible within the GL formalism where the single 
vortex will always sit in the center of the square. 
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FIG. 3. (Color online) Contour plots of the absolute value of 
the order parameter (left), the phase of the order parameter 
(center) and the LDOS at E/F = 0.57 in (c) and E/F = 
0.48 in (f) (right) for k F ^o = 2. Panels (a), (b) and (c) are for 
$/3?o = 4 and L = 1 state. Panels (d), (e) and (f) correspond 
to $/$o = 6 and L = 2 state. 

From an experimental point of view, in the absence of 
any reference energy, the zero-bias LDOS is most rel- 
evant. Here instead we will show the LDOS for the 
lowest energy vortex bound state, which for an isolated 
vortex could be found from a simple empirical formula 
E low /A = (2k F ^ a )- 1 \n{3.'S3k F ^ Q ) 0. Therefore for the 
L = 1 state we plot the LDOS at E/F = 0.57 [Fig. Etc)]. 
Note that this lowest vortex bound state is not localized 
in the vortex core but is shifted towards the center of 
the square. We attribute this to the interactions of the 
quasiparticles not only with the four deepest minima of 
the inhomogeneous order parameter but also with the 
edges and the corners - where the order parameter is also 
suppressed [see inset of Fig. Ufa)]. The effect of this in- 
teraction can also be inferred from the finite LDOS at 
the corner of the sample, next to the vortex. 

When increasing further the magnetic field, an addi- 
tional vortex enters the system (forming L = 2 state) and 
another unexpected spatial distribution is stabilized. We 
illustrate this in Figs. EJd-f) for applied flux $/$() = 6. 
The confinement seems to act strongly and vortices are 
compressed closer to each other. However, the enhance- 
ment of the order parameter in the center of the sam- 
ple due to quantum resonance prevents the two vortices 
from merging. As a consequence, vortices are squeezed 
into elliptical shapes, as a pair parallel to one of the sam- 
ple edges. This vortex configuration is as different as one 
can be from the known GL results, where the two vortices 
are always found sitting on the diagonal, or merged into 
a giant vortex, and have always an almost circular core. 
The LDOS plot [Fig. [3][f )] again reveals strong compet- 
ing interactions, different from those acting on vortices. 
For example, we see evidence of the interaction of bound 
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FIG. 4. (Color online) Contour plots of the absolute value of 
the order paramenter (a), the phase of the order parameter 
(b) and the LDOS at E/F = 0.57 in (c) for k F £o = 3 and 
$/$o = 7 (L = 3). Panel (d) shows the vortex configurations 
for the L — 4 state and their corresponding energies. 



states inside the vortex cores, since the maximum in the 
LDOS is reached between the vortices and not at the cen- 
ter of each vortex. Also, the vortex-surface interaction 
of the bound states is enhanced, leading to LDOS being 
clearly appreciable near the surface. 

To check further the influence of the length-scales 
in our sample, we also calculated the vortex states for 
kp^o = 3. In this case, the inset of Fig. [ljb) shows 
six local maxima along the diagonal, twice as many as 
in the kp£,o — 2 case. There are also strong oscilla- 
tions near the corners but relatively flatter away from 
the boundary. Again, as expected, the four-fold symme- 
try is maintained, but now there are no strong minima 
in the fourfold arrangement deep inside the sample. For 
that reason, for $/$o = 3, the L = 1 ground state is 
conventional and contains one vortex at the center of the 
square. For $/$o = 5.8, the L — 2 state is the ground 
state and although it still shows the vortex pair parallel 
to the side of the square, the shape of the vortices and 
the location of the lowest bound state LDOS peaks bring 
it closer to the conventional picture |37| . 

However, when the flux through the square is increased 
to $/$o = 7, as shown in Figs. BJa-c), the ground state 
has vorticity L — 3 and is not conventional. We in 
fact find the vortex-antivortex (v-av) molecule, similar 
to the symmetry-induced ones predicted by the GL the- 
(4 vortices with an anti- vortex in between, so 
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ory 

that the total vorticity is L = 4—1). There the size of the 
v-av molecule was found to be very small, possibly larger 
if artificial pinning centers are introduced |12l . [38j . In the 



present case, the vortex-antivortex molecule is stable over 
a wide range of fields because of the symmetry of the os- 
cillatory pattern of the order parameter due to the quan- 
tum confinement. The LDOS at E/F — 0.57, shown in 
Fig-SJc), reveals again that since the vortices are located 
closer together, the bound states are not localized in the 
cores but are extended over the whole square. 

Another found difference from earlier studies is the be- 
havior of the L = 4 state for kp^Q = 3. As shown in the 
free energy curve in Fig. 0Jd), we revealed a continu- 
ous phase transition between the configuration with the 
vortices located on the diagonal (at lower fields) and the 
configuration with the four vortices sitting near the edges 
of the square (higher fields) , a state never found within 
GL! The transition between the two four-fold symmet- 
ric states is quite peculiar and yet unseen in mesoscopic 
superconductivity - it involves the appearance of vortex- 
antivortex pairs near the center of the square [39| . As the 
field is increased, the diagonal vortices annihilate with 
the central anti- vortices, and central vortices move to the 
side location. Moreover, the L = 4 configuration with 
side vortices becomes the ground state for $/<!>o = 8.5. 
Due to the inhomogeneity of the order parameter, these 
vortices never merge into a giant vortex, contrary to the 
known GL picture for samples of smaller sizes. 

In conclusion, we found novel vortex states with un- 
conventional stability ranges and unconventional transi- 
tion between them in a superconducting square in the 
quantum limit, where significant departures from previ- 
ous works based on the GL theory are found. Experimen- 
tally, these states can be accessed through STM measure- 
ments. Additionally, we showed that competing interac- 
tions in the quantum limit for the bound states arc differ- 
ent from those for the vortices, so that the conventional 
picture of a vortex bound to lowest energy states does not 
hold. Instead we predict that the maxima in LDOS of 
the lowest energy states will be observed between vortices 
and near surfaces. These peculiar phenomena are made 
possible by strong quantum confinement, which induces 
spatial oscillations in the order parameter. Their specific 
pattern depends on the ratio of £o and Xf, which is un- 
favorable for oscillations in elementary superconductors, 
but is small enough in high-T c materials. However, to 
observe these novel states in the latter case, one should 
deal with very small samples. As an alternative, we pro- 
pose the study of a graphene flake in contact to a su- 
perconducting film, where the Fermi energy of graphene 
can be tuned by a gate. In the case of graphene on Pb, 
our calculations show that one could tune kpS,o in the 
broad range of 0.1 — 10 by shifting the Fermi energy in 
a 400x400nm flake from 0.01 to 0.1 eV above the Dirac 
point. Another system where effects of quantum confine- 
ment on vortex matter can be probed systematically are 
the optically trapped cold gases [Z(|, which are nowa- 
days extremely controllable. Further investigations will 
address the rich physics in the quantum limit, and show 
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the effects of our findings in the 3D-confined case 41 - 43| 
and multi-condensate samples but where also barri- 
ers for vortex motion across the oscillating landscape can 
be investigated for, possible use as Q-bits or other vortex 
devices [45 1 . 
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